Preprint typeset in JHEP style - HYPER VERSION 



IFUP-TH 2003/13 



A model for string-breaking in QCD 



D. Antonov* L. Del Debbio, A. Di Giacomo 

INFN-Sezione di Pisa, Universitd degli studi di Pisa, Dipartimento di Fisica "E. 
Fermi", Via Buonarroti, 2 - Ed. B-C - 56127 Pisa, Italy 
E-mail: 



antonov, Idd, digiaco@df.unipi.it 



Abstract: We present a model for string breaking based on the existence of chromoelectric 
flux tubes. We predict the form of the long-range potential, and obtain an estimate of the 
string breaking length. A prediction is also obtained for the behaviour with temperature 
of the string breaking length near the deconfinement phase transition. We plan to use this 
model as a guide for a program of study of string breaking on the lattice. 
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I. Introduction 

In pure gauge theories, the potential energy of a static qq pair grows with the quark- 
antiquark distance r as 

V{r) = (jr (1.1) 

where a is known as string tension. From Regge phenomenology, a ~ (440 MeV)^. Equa- 
tion ( |1.1| ) is confirmed by numerical simulations on the lattice, see e.g. [|l|], or and 
references therein for recent calculations. 

In full QCD, due to the presence of dynamical quarks, string breaking is expected to 
occur. Equation ( p!j| ) yields an accurate description of the potential up to some distance, 
above which the system prefers to convert the energy into qq pairs, or mesons. String 
breaking has been observed on the lattice, see e.g. ^ ^, and its behaviour at finite 
temperature around the deconfining phase transition has been studied numerically Q. 

In this paper, we construct a dynamical model of string breaking. Lattice simulations 
have shown that the chromoelectric field in the presence of a static quark-antiquark pair 
looks like an Abrikosov flux tube Q of transverse size r_L, with r_L ~ 0.2 0.3 fm g, |lO|, 

II, O, O, 111], as conjectured since a long time |15]. This picture is consistent with a dual 



superconductor mechanism of confinement |16, 17|, the superconductor being at the border 
between type-I and type-II. We shall interpret string breaking as due to the production of 
light quark-antiquark pairs by the chromoelectric field of the dual Abrikosov-Nielsen-Olesen 
flux tube. We plan to verify this idea by an extended program of lattice simulations. In this 
paper, we discuss the general features of the model and present some analytic calculations 
in limiting cases. 

A similar phenomenon occurs in electrodynamics when one considers the potential 
inside a condenser. When the distance r between the plates is increased, at constant 
electric field, the energy increases linearly with r, up to some distance where the system 
becomes unstable because of e"'"e~ pair production. 
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For a constant static electric field the rate of production of spin-1/2 particle pairs per 
unit volume and unit time is given by [|l^ ]: 



^ ' n=l 

In QCD the analog of the condenser is the flux tube, electron-positron pairs are replaced 
by quark-antiquark pairs, the coupling constant e by the strong coupling constant g, and a 
factor NcNf must be added to properly take into account the number of degrees of freedom. 
Hence, 

^QCD = 2N^Nf^^ ^2 ' (1-3) 

^ ' n=l 

where if is a constant static chromoelectric field. However, free quarks do not exist in QCD, 
and therefore m in this formula is some mass parameter to be determined, which can be 
either the constituent quark mass, or a typical hadron mass (e.g., rriTr). On the lattice, 
quark masses can be changed at will, and m can be understood as a function of them. 

In this work, in order to understand the relevant physics, we shall study the simplified 
case, in which the chromoelectric field is averaged over the transverse direction and then 
considered as uniform in space. A partial physical justification is that, due to the higher 
modes of the string, the tube configuration will fiuctuate with a frequency that is much 
higher than the rate of pair production. More precise calculations taking into account 
the transverse size of the flux tube and the space dependence of the electric fleld can 
be performed and will be presented elsewhere. The use of Eq. ( |1.2[ ), which relies on the 
assumption that the field is constant in space, will in any case be a good approximation if 
the transverse size of the tube is larger than 1/m. The volume that we consider is the 
physical three-dimensional volume of the flux tube, of cross section S, multiplied by a time 
equal to the length of the tube r in natural units. This is a natural choice. Indeed, in order 
to measure the force between the heavy qq pair, we need at least the time to propagate the 
interaction between the two particles. Also this choice can be tested on the lattice. The 
interaction can be described by an elastic potential if no pair has been produced in the 
meantime. Hence, the long-range part of the potential is given by: 

V{r) =are"''''^'" (1.4) 

i.e. the linear potential ar multiplied by the probability that no pairs are created in the 
volume rj. Clearly, V{r) has a maximum at: 

(1.5) 



V2Sw 



which yields an indication of the string breaking distance. Models for hadronic processes 
based on similar ideas exist in the literature (see e.g. p9[| ). 
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2. Explicit models 



We shall assume for the average chromoelectric field E the rms value of the chromoelectric 
field on the cross section of the tube: 



At the border between type-I and type-II supeconductor (the so-called Bogomol'nyi limit 
||20|]), {S"^) can be computed explicitely. The flux tube is the solution of the dual 
Landau-Ginzburg equation, obtained from the Lagrangian: 

C = -\f' + \D^\^-V{^). (2.1) 

Here F^^ = df,B^ - dyB^, D^^p = {d^ - igmB^)ip, V{ip) = | [\^p\'^ - v^)^, is the dual 
gauge field, gm is the magnetic charge and, in the broken phase, v ^ 0. The Bogomol'nyi 
limit corresponds to mu = my = V^gm v- In this limit, the solutions to the equations of 
motion are known and the string tension, defined as the integral of the energy density 



H over transverse directions, can be computed exactly |20|: 

a = [ d'zn{z) = ^. (2.2) 

J dm 

Using the Dirac quantization condition, ggm = 27r, Eq. (|2.2| ) becomes ar'j_ = Og. 

Using the solution of Ref. [^] for a tube carrying one unit of flux, we compute the 
average of the electric field squared: 

{£') = C^o,^, (2.3) 
vr 

where Ceog — 0.18. The exponent in Eq. ( |l.3D becomes then 

VO.lSggmO- 

and if ggm = 27r, as expected in a dual superconductor, then y ~ 0.7. The sum of the 
terms of order higher than one in Eq. (|1.3| ) can be estimated to be ~ 2/^/4, that for the 
worst case of m = allows one to keep only the first term in the sum with an accuracy 
of ~ 17%. The coefficient w can finally be written as: 



iV.iV;^ ^ae-'^"^'/2^^ (2.5) 



w 



^2 S 



Note that the potential in this analysis is determined up to a constant term, depending 
on the physics at shorter distances. Furthermore, V{r') only depends on the string tension 
a and the mass m. The potential obtained by inserting Eq. ( p.5| ) into Eq. ([1.4]) can be 
checked on the lattice; the dependence on rr? can be tested by comparing the potential 
while varying the bare quark mass. 
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Equation ( |l.5D yields the following expression for the position of the maximum: 

vr 



^OMN.Nfasay 



(2.6) 



Note that since asNc is constant at large Nc, r is A'^c-independent in this limit. Next, 
using Q!s(l GeV) ~ 0.5, Nc = 3, Nf = 2, one obtains from Eq. (2.6) the following estimates 



f{m = 200 MeV) ~ 2.0 fm, r{m = 300 MeV) ~ 3.3 fm. There is a degree of arbitrariness 
in the choice of the scale at which ag should be computed, so that the latter has to be 
considered as a parameter in Eq. (p. 61) . In practice § ^, ^ , string breaking is detected 
on the lattice as a deviation of the gg-potential from the straight line V = ar, bigger than 
the numerical error. Let e be the value of such a relative error. From Eq. (|1.4[) it follows 



1 _ e-^S"- ^ g 



or the breaking distance r^: 



For instance, for e ~ 20%, ~ 0.28f. 

The existence of the maximum of the potential and the dependence of its position on 
m can also be tested directly in numerical simulations. In all this analysis the flux tube is 
assumed to be that of the quenched theory. Production of gluon pairs is neglected, because 
the threshold would be given by the lightest glueball mass, which is much larger than the 
typical hadron mass. This is also the reason why at the one-loop order the formula ( |1.2| ) 
can be translated to QCD in the way this has been done in Eq. ( |1.3| ). 

The same analysis can be performed in the extreme type-II regime (else called the 
London limit), Iuk ^ 1, where k = mH/mv is the Ginzburg-Landau parameter with 
niH = \/2At;, my = \/2gmV. The only difference of this limit from the Bogomol'nyi one is 
that the coefficient l/\/0.18 ~ 2.36 in Eq. ( |2.4| ) should be replaced by I.ISIuk. Indeed, in 
the London limit Eq. ( |2.3D is replaced by 

{S^) = (2.7) 
m K vr 

where CLond — 0.71. For reasonable values of Iuk (for instance, Iuk ~ 4), there is no 
significant change with respect to the Bogomol'nyi case. 



3. Temperature dependence 

More information can be extracted from Eq. ( |1.5D if one considers the long-range potential 
at finite temperature, as the deconfinement temperature is approached. 

If the mass in Eq. (|1.2| ) is the constituent quark mass, it stays finite at the transition, 
while the string tension vanishes with a critical exponent: 

ar.{l-T/T,Y (3.1) 
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where an effective exponent = 1/3 can be used for a first order pfiase transition. In this 
case, f behaves as ^"""^"^ _ The approximation in Eq. (^]^) is improved as we approach the 
phase transition. 

If instead the mass in Eq. ( |1.2| ) is the pion mass, it vanishes at the chiral point as: 

m2 ~ (1 - T/T,)\ (3.2) 
where the exponent 7 is determined by the universahty class of the transition. For instance. 



if the symmetry breaking is described by an 0(4) RG fixed point |22|, 7 ~ 1.44. In this case, 
y approaches one at the critical temperature, and the string breaking distance grows as 
\l\fo as the temperature is increased. One should also remark that Eq. (|2.2| ) implies that 
the transverse size of the flux tube r_L is inversely proportional to \fa. When approaching 
Tc, r\ goes large as (1 — T/Tc)~'^ ■ Since v = 1/3, r\m?' ~ (1 — T/TcY'^ goes to zero, 
implying that the Schwinger formula can no longer be applied when T is very near T^. 
These two scenarios can be checked against data from lattice simulations. 

A qualitative comparison with the data of Ref. Q for the string breaking length at finite 
temperature is roughly consistent with the behaviour predicted by combining Eq. ( |2.6D with 



Eqs. (|3.1|) and (3.2). A more quantitative analysis requires more precise and systematic 



numerical simulations. 



4. Concluding remarks 

We have developed a model for string breaking based on the existence of chromoelectric 
flux tubes. The form of the long-range potential is predicted, which leads to an estimate 
of the string breaking length; the relevant parameters are the string tension and the mass 
m, which is to be determined phenomenologically. A prediction is also obtained for the 
behaviour of the string breaking length near the deconfinement phase transition. Our 
predictions will be tested against the results of lattice simulations. In particular, the quark 
masses in lattice simulations can be varied, and m can be understood as a function of them. 
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